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FOREWORD 


This  user's  manual  documents  the  Navy  Energy  Siting  (NES)  code 
developed  by  Acurex  Corporation,  Energy  and  Environmental  Division, 

Mountain  View,  California  for  the  Civil  Engineering  Laboratory,  Naval 
Construction  Battalion  Center,  Port  Hueneme,  California.  This  work  was 
performed  on  Contract  N68305-78-C-0009. 

This  user's  manual  represents  Volume  II  of  the  Final  Report.  It 
provides  the  information  necessary  to  understand  input  requirements  and 
program  macro  logic  of  the  NES  code. 
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SECTION  1 


INTRODUCTION 


1.1  PURPOSE 

The  Navy  Energy  Siting  (NES)  code  was  developed  to  determine,  at 
various  naval  sites,  the  optimum  mix  of  alternate  energy  and  commercially 
purchased  energy  that  will  meet  all  the  site's  energy  demands.  The 
optimum  mix  meets  the  energy  demands  within  given  site  constraints  at  the 
lowest  cost. 

The  code  handles  three  types  of  energy  demands:  electricity,  space 
heating  and  hot  water  (combined),  and  process  steam.  Alternate  energy 
models  are  designed  to  satisfy  a  single  demand,  with  the  exception  of  a 
cogeneration  model  which  can  satisfy  two  energy  demands  concurrently. 

1.2  PROCESSING  PERFORMED  8Y  THE  OPTIMIZATION  PROGRAM 

The  gode  uses  a  gradient  projection  method  to  determine  the  optimum 
mix  of  systems  that  minimize  energy  cost.  The  gradient  projection 
routines  were  provided  by  the  Civil  Engineering  Laboratory  while  the 
alternate  energy  models  and  supporting  routines  were  developed  by  Acurex. 

During  each  iteration,  the  optimization  routine  selects  a  mix  of 
alternate  energy  systems  which  reduces  total  energy  costs  relative  to  the 
system  mix  selected  during  the  previous  iteration.  Total  energy  cost  is 
determined  by  totaling  the  annualized  life-cycle  costs  of  each  alternate 
energy  system  and  commercial  energy  costs.  The  amount  of  commercial 


energy  that  must  be  purchased  (electricity,  space  heating,  hot  water,  and 
process  steam)  is  calculated  by  deducting  the  energy  produced  by  the 
alternate  energy  sources  from  the  total  demand.  This  iterative  procedure 
continues  until  a  minimum  annual  cost  is  determined. 

1.3  RESTRICTIONS  AND  LIMITATIONS 

Restrictions  and  limitations  associated  with  use  of  the  NES  code 
are  listed  below: 

•  The  user  is  restricted  to  three  energy  demands  and  the 
alternate  energy  models  currently  programmed  in  the 
optimization  code.  The  list  of  alternate  energy  models  is 
given  in  Table  1-1. 

•  When  selecting  starting  values  for  the  alternate  energy  models 
the  user  should  pick  values  which  do  not  exceed  energy  demands, 
site  area,  and  coal  and  refuse  constraints  on  the  first 
Iteration.  If  the  demands  are  greatly  exceeded  the 
optimization  process  does  not  converge. 

•  The  user  should  select  reasonable  upper  and  lower  bounds  on  the 
number  of  systems  for  each  alternate  energy  model  because  the 
difference  between  the  bounds  determines  the  step  size  for  the 
first  iteration. 
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TABLE  1-1.  ALTERNATE  ENERGY  MODELS  LIST 


Subroutine 

Name 

Model 

Name 

Model 

No. 

Description 

EVAL01 

SLTHTG 

1 

Solar  thermal  for  space  heating  and 
hot  water. 

EVAL02 

RDFHTG 

2 

Refuse  derived  fuel  for  space  heating  and 
hot  water. 

EVAL03 

RDFSTM 

3 

Refuse  derived  fuel  for  process  steam. 

EVAL04 

RDFELE 

4 

Refuse  derived  fuel  for  electricity. 

EVAL05 

FBCHTG 

5 

Fluidized  bed  combustion  for  space 
heating  and  hot  water. 

EVAL06 

FBCSTM 

6 

Fluidized  bed  combustion  for  process 
steam. 

EVAL07 

FBCELE 

7 

Fluidized  bed  combustion  for 
electricity. 

EVAL08 

GEO STM 

8 

Geothermal  for  process  steam. 

EVAL09 

GEOELE 

9 

Geothermal  for  electricity. 

EVALIO 

WD5 

10 

5  kW  wind  generator  for  electricity. 

EVAL11 

WD200 

11 

200  kW  wind  generator  for  electricity. 

EVAL12 

WD1500 

12 

1500  kW  wind  generator  for  electricity. 

EVAL13 

PHVELE 

13 

Photovoltaic  for  electricity. 

EVAL14 

CCLHTG 

14 

Conventional  coal  for  space  heating  and 
hot  water. 

EVAL15 

CCLSTM 

15 

Conventional  coal  for  process  steam. 

EVAL16 

CCLELE 

16 

Conventional  coal  for  electricity. 

EVAL17 

CCLCOG 

17 

Conventional  coal  for  cogeneration  of 
process  steam  and  electricity. 
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SECTION  2 


HOW  TO  USE  THE  CODE 


The  input  necessary  to  run  the  NES  code  is  broken  into  three 
categories:  site  data,  model  data,  and  run  specification  data.  The  deck 
structure  is  depicted  in  Figure  2-1. 

The  site  input  data  is  site  dependent  data  which  changes  from  site 
to  site.  This  data  includes  site  energy  demand  profiles,  site  insolation 
profiles,  site  wind  profiles,  and  general  site  inputs  such  as  commercial 
energy  costs,  conmercial  energy  purchase  limits,  and  availability  of  coal 
and  refuse. 

The  model  input  data  is  model  dependent  data  required  by  each 
alternate  energy  model  to  compute  energy  produced,  cost,  and  area 
required.  This  data  includes  such  parameters  as  capital  cost,  fuel  cost, 
efficiency,  and  area  factors.  A  complete  list  of  parameters  for  each 
alternate  energy  model  is  given  in  Volume  I,  Appendix  A.  Actual  data 
required  is  described  subsequently  in  Tables  2-6  to  2-9. 

The  run  specification  input  data  specifies  which  alternate  energy 
models  are  to  be  considered  in  a  particular  run. 

2.1  SITE  DATA 

For  simplicity,  the  site  input  data  consist  of  a  fixed  number  on 
input  records  (cards).  The  code  assumes  that  all  site  data  is  present. 
Therefore  all  data  cards  must  be  present,  although  some  may  be  blank  if 
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the  data  is  not  to  be  used  in  a  particular  run.  For  example,  if  the  user 
is  considering  only  electrical  demand,  then  data  for  the  other  two  demands 
need  not  be  input,  but  the  blank  data  cards  must  be  present. 

The  site  data  is  divided  into  seven  subsets  and  the  order  in  which 
they  occur  in  the  code  is  as  follows: 

•  PROBLEM  TITLE  SET  (3  cards) 

•  ENERGY  DEMAND  SETS  (36  cards  each) 

•  COMMERCIAL  ENERGY  SET  (4  cards) 

•  GENERAL  DATA  SET  (6  cards) 

•  INSOLATION  DATA  SETS  (36  cards  each) 

•  WIND  VELOCITY  DATA  SET  (36  cards) 

•  GEOTHERMAL  DATA  SET  (1  card) 

Each  subset  is  described  in  detail  below. 

Unless  otherwise  specified,  the  user  may  assume  that  all  numeric 
input  data  is  contained  in  floating  point  fields  of  width  10  columns,  in 
which  E  format  is  acceptable  (E10.0).  Furthermore,  it  is  assumed  that  all 
yearly  input  profiles  begin  in  January  and  end  in  December. 

2.1.1  Problem  Title  Set 

This  data  set  consists  of  three  cards  describing  the  case.  Columns 
1  through  78  are  available  on  each  card  (Figure  2-2)  and  the  information 
is  printed  at  the  beginning  ef  the  run.  Format  for  these  cards  is  (13A6). 

2.1.2  Energy  Demand  Sets 

Alternate  energy  models  compete  in  three  energy  demand  sectors: 
electrical,  space  heating  and  hot  water  (combined),  and  process  steam. 

Each  energy  demand  set  consists  of  36  cards  which  contain  12 
monthly-average  days  of  hourly  demand  data.  Each  card  represents  8  hours 
of  demand.  Therefore  three  cards  make  up  a  daily  profile  of  24  hours 
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cial  energy  card  sets 


(Figure  2-2).  All  cards  must  be  present  for  all  three  demands  but  only 
the  demand  sets  being  modeled  need  to  have  actual  data  input. 


Format  for  each  card  in  the  energy  demand  sets  is  (8E10.0). 

Table  2-1  indicates  the  sequence  and  units  of  the  demand  inputs. 

2.1.3  Commercial  Energy  Set 

This  set  of  four  cards  (only  the  first  three  cards  are  used) 
contains  information  about  the  commercial  energy  purchased  to  satisfy  the 
three  energy  demands.  The  cards  are  input  in  the  same  order  described  for 
the  energy  demand.  The  fourth  is  left  blank.  Each  card  is  identically 
formatted  and  described  in  Figure  2-2  and  Table  2-2. 

2.1.4  General  Data  Set 

This  data  consists  of  six  cards  containing  general  information 
about  the  site.  A  summary  of  this  information  and  its  use  follows.  Card 
formats  are  described  in  Table  2-3  and  in  Figure  2-3. 

Card  1 

The  quantity  of  coal  (tons/day)  available  for  the  coal  combustion 
energy  models  is  limited.  Also,  capital  cost  of  coal  combustion  systems 
depends  upon  the  sulfur  content  of  the  coal.  The  amounts  (tons/day)  and 
the  particular  type  of  coal  (high  or  low  sulfur)  available  at  a  site  is 
input  on  card  1. 


TABLE  2-1.  ENERGY  DEMAND  SETS 


Number 

Demand 

Units  of  Input 

1 

Electrical 

MWh's 

2 

Space  heating  and  hot  water 

MBtu's 

3 

Process  steam 

MBtu's 
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TABLE  2-2.  COMMERCIAL  ENERGY  SET3 


Column 

Input 

Units 

1  -  10 

Commercial  energy  cost 

($/MWh)  Electrical 
($/MBtu)  Space  heating 
($/MBtu)  Process  steam 

11  -  20 

Differential  inflation  rate 

(%)  expressed  as  a  fraction 

21  -  30 

Maintenance  and  operation  cost 

(it)  expressed  as  a  fraction 

31  -  40 

Efficiency 

(X)  expressed  as  a  fraction 

41  -  50 

Capital  cost 

(S) 

51  -  60 

Purchase  limit 

(MWh/hr)  Electrical 
(MBtu/hr)  Space  heating 
(MBtu/hr)  Process  steam 

61  -  70 

_ 

Purchase  limit 

(MWh/yr)  Electrical 
(MBtu/yr)  Space  heating 
(MBtu/yr)  Process  steam 

aFormat  (7E10.4) 


TABLE  2-3.  CARD  FORMAT  FOR  GENERAL  DATA  SET 


Card/Column 

Format 

Input 

Card  1 

11  -  20 

E10.0 

Coal  auantity  (tons /day) 

26  -  32 

Card  2 

A6 

Coal  quality  ("high"  or  "low") 

11  -  20 

Card  3 

E10.0 

Land  area  available  (ft^) 

11  -  20 

Card  4 

E10.0 

Refuse  available  (tons/day) 

1  -  20 

Card  5 

E10.0 

Discount  interest  rate 
(%  expressed  as  a  fraction) 

1  -  80 

Card  6 

8E10.0 

Ambient  temperature  (°F) 
January  -  August 

1  -  40 

4E10.0 

Ambient  temperature  (°F) 
September  -  December 

2-7 


inn 

Him 

nm 

Inn 

■I  IH 
!■■■ 


IB 


2-8 


Card  2 


2 

Card  2  contains  the  land  area  (ft  )  available  for  siting 
alternate  energy  systems. 

Card  3 

Similar  to  coal,  refuse  is  limited  at  each  site.  The  amount 
available  (tons/day)  is  input  on  c<nd  3. 

Card  4 

The  discount  (interest)  rate  (percent  expressed  as  a  fraction)  used 
in  the  uniform  annual  cost  function  (see  Volume  I,  Section  2.3)  is  given 
on  card  4. 

Cards  5  to  6 

These  cards  contain  the  monthly  average  temperatures  (°F) 
representing  one  year  of  data.  These  data  are  used  by  Subroutine  FCHART 
to  determine  the  amount  of  energy  supplied  by  the  solar  thermal  model. 

2.1.5  Insolation  Data  Sets 

There  are  two  insolation  data  sets  (36  cards  each):  one  for  the 
solar  thermal  model  and  one  for  the  photovoltaic  model.  Each  data  set 
contains  hourly  insolation  profiles  for  12  average  days  representing 
12  months  of  the  year.  Identical  to  demand  data  card  (see  Figure  2-2), 

each  card  contains  8  hours  of  data  with  three  cards  comprising  a  24-hour 

? 

profile.  The  units  of  input  are  (Btu/ft  /hr).  The  format  for  each  card 
is  (8E10.0). 

2.1.6  Wind  Data  Set 

The  wind  data  set  consists  of  36  cards  representing  the  hourly  wind 
velocity  profiles  for  the  site.  The  set  contains  profiles  for  12  monthly 
average  days,  24  hours  for  each  day  (Figure  2-2).  Each  card  contains 
8  hours  of  data,  formatted  (8E10.0).  Units  of  input  are  (MPH). 
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2.1.7  Geothermal  Data  Set 

Three  inputs  are  needed  for  the  geothermal  model:  geothermal  pool 
quality  (LIQUID  or  VAPOR),  geothermal  pool  size  (MBTU),  and  geothermal 
pool  temperature  (DEG  C).  The  format  for  this  data  is  given  in  Table  2-4 
and  illustrated  in  Figure  2-4. 

2.2  MODEL  DATA  SETS 

The  model  data  sets  consist  of  input  for  each  alternate  energy 
model.  Unlike  the  site  data,  model  data  should  be  input  only  if  the  model 
is  a  candidate  model  in  the  run.  The  order  of  input  determines  the 
sequence  in  which  the  models  are  called  to  satisfy  demands.  The  alternate 
energy  models  currently  available  in  the  code  were  listed  in  Table  1-1. 

The  data  set  required  by  each  model  consist  of  five  cards.  The 
format  for  these  data  sets  was  standardized  as  much  as  possible  to 
simplify  input  procedures.  The  standard  format  is  illustrated  in 
Figure  2-5. 

The  user  should  review  the  alternate  energy  models  described  in 
Volume  I,  Appendix  A  before  using  the  tables.  The  first  card  of  each 

TABLE  2-4.  CARD  FORMAT  FOR  GEOTHERMAL  MODEL 


Column 

Format 

Input 

1  -  6 

A6 

"LIQUID  or  VAPOR" 

Vapor  sets  geothermal  efficiency  at  20% 

Liquid  sets  geothermal  efficiency  based  on  pool 
temperature 

Efficiency  =  2.308E-04  *  Temp  (°C)  +  0.0392 

11  -  20 

E10.0 

Pool  size  (MBtu) 

21  -  30 

E10.0 

Pool  temperature  (°C) 
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gure  2-5.  Card  format  for  the  model  data  set. 


model  data  set  contains  the  abbreviated  name  for  the  energy  model.  The 
abbreviated  names  were  given  in  Table  1-1.  The  remaining  four  cards  in 
the  model  data  set  are  described  in  Tables  2-5  through  2-8.  The  tables 
are  set  up  in  a  matrix  format  where  each  table  describes  a  particular  card 
in  the  data  set.  The  rows  of  the  matrix  describe  a  particular  field  on 
the  card.  The  columns  in  the  matrix  represent  each  of  the  alternate 
energy  models. 

These  columns  are  labeled  by  the  abbreviated  names  of  the  models. 
The  body  of  the  matrix  contains  the  units  and/or  description  of  inputs  for 
that  model.  All  fields  described  in  these  tables  are  formatted  for  input 
(E10.0).  For  example,  for  card  2  of  SLTHTG  (solar  heating  model),  four 
parameters  should  be  input:  two  capital  cost  factors,  a  maintenance  and 
operating  cost,  and  an  exponent. 

The  code  recognizes  the  end  of  the  alternate  energy  model  data  set 
when  it  sees  ENDMDL  (end  model)  in  columns  1-6.  This  causes  the  code  to 
stop  searching  for  model  data  and  to  look  for  a  run  specification  data  set 
(Figure  2-1). 

2.3  THE  RUN  SPECIFICATION  DATA  SET 

The  run  specification  data  set  identifies  the  conditions  of  the 
run:  the  energy  demands  and  the  alternate  energy  models  to  be 
considered.  This  data  set  consists  of  four  cards.  Card  1  specifies  the 
demands  to  be  considered,  while  cards  2  and  3  specify  models  to  be  used. 
Card  4  specifies  a  print  option. 

The  user  must  provide  all  necessary  input  for  the  energy  demands 
and  alternate  energy  models  being  considered.  If  an  alternate  energy 
model  is  requested  and  no  data  was  input,  the  program  will  terminate  with 
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TABLE  2-5.  DESCRIPTION  OF  REQUIRED  DATA  AND  CARD  FORMAT 
FOR  MODEL  COST  DATA 
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a  diagnostic  error  message.  Similarly,  if  a  demand  or  alternate  energy 
model  name  is  not  recognized,  the  program  will  terminate  with  a  message. 
2.3.1  Card  Formats 
Card  1 

Card  1  specifies  the  energy  demand  to  be  considered.  The  energy 
demand  six  character  names  appear  left  justified  in  any  of  three  fields: 
columns  1-6,  columns  8  -  13,  columns  15  -  19. 

The  energy  demand  names  are  as  follows: 

•  ELECTR  Electrical 

•  SPCHTG  Space  heating  and  hot  water 

•  PROSTM  Process  steam 

Cards  2  and  3 

The  alternate  energy  model  names  appear  left  justified  in  any  of  10 
fields  on  each  card.  The  order  of  model  names  input  determines  the 
calling  sequence  of  the  models  considered.  The  fields  begin  in  columns  1, 
8,  15,  22,  29,  36,  43,  50,  57,  64.  Only  those  alternate  energy  models 
being  considered  appear  on  these  cards.  The  six  character  model  names  are 
listed  in  Table  1-1.  Both  cards  must  be  input  even  if  one  is  blank. 

Card  4 

Print  option  flag  is  specified  on  this  card. 

Table  2-9  describes  the  output  generated  for  various  values  of  the 
print  option  flag.  A  print  option  flag  of  a  particular  value  will 
generate  information  for  that  value  as  well  as  information  for  flags  of 
lower  value.  For  example,  if  the  user  specifies  Flag  1,  he  will  receive 
all  the  output  associated  with  the  Flags  -1,  0,  and  1.  Needless  to  say, 
considerable  output  is  produced  for  a  print  option  flag  other  than  -1. 
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TABLE  2-9.  PRINT  OPTION  FLAG  OUTPUTS 


Flag 

Printed  Output 

-1 

Input  data. 

Initialization  data. 

Results  summary. 

0 

Overall  status  of  the  optimization  process. 

Status  after  each  evaluation  of  alternate  energy  models. 

Partial  derivatives. 

Objective  function. 

1 

Point  to  be  projected  onto  constraints. 

2 

Status  of  energy  demands  after  evaluation  of  alternate  energy 
mode  1 s . 

Values  of  various  constraints. 

Matrixes,  vectors  and  solutions  used  during  optimization. 

Column  Format  Input 

1-2  12  PRINT  OPTION  FLAG 

(-1,  0,  1  or  2) 
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SECTION  3 


DESCRIPTION  OF  OUTPUT 

The  printed  output  generated  by  the  Navy  Energy  Siting  code  fall 
into  the  following  categories:  input,  summary  of  results,  debug  and 
diagnostics. 

All  site  model  and  run  specification  input  data  is  printed  at  the 
beginning  of  the  run.  A  summary  of  optimum  results  is  printed  at  the  end 
of  the  run.  As  an  option,  a  debug  printout  monitors  the  status  of  the 
optimization  process  (Flag^O).  This  debug  printout  indicates  the  status 
of  each  alternate  energy  model,  costs  of  alternate  energy  produced, 
cost  of  commercial  energy  purchased,  and  total  cost  of  energy  supplied 
(objective  function).  The  debug  option  generates  considerable  output  and 
is  not  recommended  unless  it  is  necessary  to  monitor  the  optimization 
process . 

In  addition,  a  set  of  diagnostic  error  messages  were  incorporated 

into  the  code.  These  are  listed  below. 

3.1  DIAGNOSTIC  MESSAGES 

"***ERR0R  DEATH***CARD  READ  ERROR" 

Usually  a  result  of  improper  deck  set  up,  a  card  missing, 
or  a  card  out  of  sequence. 

"***ERR0R  DEATH***PREMATURE  END  OF  DATA" 

End  of  deck  reached  before  all  data  is  read.  Check  deck 
for  missing  data. 
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"***ERROR  DEATH***MODEL  NOT  RECOGNIZED" 

When  trying  to  input  model  data,  the  model  name  was  not 
recognized. 

"***ERR0R  DEATH***DATA  FOR  A  MODEL  INPUT  TWICE" 

Self-expl anatory. 

“NO  DATA  FOR  MODEL" 

A  model  has  been  requested  for  which  no  data  has  been  input. 
"(Name)  IS  NOT  IN  THE  LIST" 

Check  run  specification  data.  A  demand  or  model  requested 
is  not  recognized. 

"NO  MATCH  ON  THE  n  VALUES" 

Check  run  specification  data,  no  demands  or  models 
recognized. 

"ERROR  ON  INPUT" 

Check  run  specification  data. 
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SECTION  4 


PROGRAM  LOGIC  DESCRIPTION 


This  section  presents  a  description  of  the  Navy  Energy  Siting  (NES) 
code  structure  and  logic  flow.  Included  in  this  section  is  a  flow  chart 
describing  the  major  processes  performed  by  the  code. 

The  Navy  Energy  Siting  code  is  driven  by  the  main  program  ENSITE, 
which  calls  various  subroutines  designed  to  handle  the  major  processes 
performed.  The  code  can  be  broken  into  three  major  components: 
input/output,  optimization  and  the  set  of  alternate  energy  models. 

The  input  process  consists  of  reading  three  types  of  data;  site 
data,  model  data  and  the  run  specification  data.  All  input/output 
variables  are  stored  in  labeled  common  blocks,  which  allow  subroutines 
easy  access  to  large  blocks  of  information. 

The  optimization  portion  of  the  program  consists  of  a  set  of 
subroutines  supplied  by  the  Civil  Engineering  Laboratory.  These  routines 
use  a  gradient  projection  technique  to  determine  the  optimum  mix  of 
alternate  energy  systems. 

The  alternate  energy  models  are  the  individual  subroutines  designed 
to  model  the  various  alternate  energy  systems. 

4.1  INPUT/OUTPUT 

The  main  program  ENSITE,  as  shown  in  Figure  4-1,  calls  three 
subroutines  to  load  the  input  data:  INSITE,  INMODL,  RUNSPC.  The  data  are 


ENSITE 


START 


INSITE 


INPUT  SITE  DATA 


INHODL 


INPUT  MODEL  DATA 


RUNSPC 


INPUT  RUN 
SPECIFICATION 


OPTIMIZATION 


OPTMIZ 


RAPUP 


ORDER  MODELS 
DETERMINE  EXCESS 
ENERGY 


PRINT  SUMMARY 


PRTSUM 


UPDATE 


USE  MAXIMUM  NUMBER  OF 
SYSTEMS  AND  PRINT  SUMMARY 
FOR  COST  COMPARISONS 


Figure  4-1.  Navy  Energy  Siting  code. 


stored  in  labeled  common  blocks  accessible  to  all  routines.  Subroutine 
INSITE  reads  the  site  dependent  input  data  and  stores  it  in  the  labeled 
common  SITE  and  FUELSC.  The  site  data  consists  of  the  demand  data: 
electricity,  heating,  process  steam,  insolation  and  wind  velocity 
profiles.  Also  included  are  geothermal  reservoir1  data,  commercial  energy 
costs,  inflation  rates,  and  crunch  factors. 

Subroutine  INMODL  reads  the  model  input  data  and  stores  the  input 
in  the  labeled  common  block  MODEL.  The  model  input  consists  of  data 
necessary  for  each  alternate  energy  model.  This  data  includes  capital 
cost  factors,  fuel  costs,  fuel  inflation  rates,  maintenance  and  operating 
costs,  exponents,  transportation  cost,  revenue  from  recovered  material 
(RDF),  efficiency  factors,  load  factors,  area  factors,  minimum  and  maximum 
number  of  systems,  and  starting  values  for  number  of  systems.  Not  every 
alternate  energy  system  requires  all  the  input  parameters  mentioned,  but 
all  parameters  are  available  to  each  model. 

Suo'outine  RUNSPC  reads  input  data,  specifying  which  energy  demand 
is  to  be  satisifed  and  which  alternate  energy  models  are  to  be  considered. 

The  output  processes  are  handled  by  a  wide  range  of  subroutines, 
each  designed  to  print  a  variety  of  information.  All  site,  model  and  run 
specification  data  are  printed.  During  the  iterative  steps  of  the 
optimization  process,  various  levels  of  information  are  also  printed 
indicating  the  progress  of  the  optimization  program.  Once  the 
optimization  is  complete,  there  is  a  summary  printout  that  lists  the  mix 
of  systems  used,  costs  and  energy  produced. 

4.2  OPTIMIZATION 

The  optimization  process  uses  the  gradient  projection  technique  to 
minimize  total  energy  costs.  Subroutine  OPTMIZ  (NONLIN)  initiates  the 
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program  by  calling  subroutine  INPUT  which  loads  some  fixed  inputs  and  the 
minimum,  maximum  and  starting  points  (number  of  systems)  for  each  model 
(see  Figure  4-2).  This  information  is  stored  in  labeled  common  LINCOM. 

The  optimization  process  calls  EVAL  to  compute  alternate  energy  cost  and 
commercial  energy  cost. 

Subroutine  CSUM  and  PARTL  are  also  called  during  the  optimization 
process  to  check  constraints  and  compute  the  partial  derivatives  of  the 
constraints.  Once  the  costs  and  value  of  constraints  are  determined,  the 
optimization  routines  adjust  the  mix  of  alternate  energy  systems 
accordingly,  to  minimize  costs  and  not  violate  the  various  constraints. 

This  process  is  repeated  until  a  minimum  is  reached. 

4.3  EVAL  (ALTERNATE  ENERGY  SYSTEMS) 

The  EVAL  Subroutine  calls  the  individual  alternate  energy  models 
(EVALn)  and  the  commercial  energy  routine  (COSTCE)  to  determine  total 
cost  to  satisfy  all  energy  demands  (Figure  4-3).  Each  alternate  energy 
model  (EVALn),  given  the  number  of  systems  to  evaluate  (X^,  computes 
energy  produced,  area  required  for  systems,  and  cost  of  energy  produced. 

This  information  is  passed  back  to  the  main  subroutine  EVAL. 

To  determine  the  demand  for  commercial  energy,  EVAL  calls  the  LOADDM 
subroutine  before  calling  the  alternate  energy  models.  This  moves  the  energy 
demands  stored  in  common  to  the  working  array  DMANDW.  Each  model  called 
produces  energy  on  an  hourly  basis,  decrements  the  working  demand  array 
(energy  produced  is  subtracted  from  the  demand);  the  remaining  energy  demand 
is  the  commercial  energy  requirement.  EVAL  then  calls  the  COSTCE  subroutine 
to  determine  the  cost  of  purchasing  conmercial  energy  that  meets  the  demand 
not  satisfied  by  the  alternate  energy  models.  Finally,  EVAL  returns  the  total 
cost  to  satisfy  the  annual  energy  demand  to  the  OPTMIZ  driver. 
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Figure  4-2.  OPTMIZ  flow  chart. 
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4.4  COMMON  BLOCK  STORAGE 


Most  of  the  communication  between  subroutines  is  through  labeled 
common  blocks.  In  this  section,  the  nature  of  the  data  contained  in  each 
labeled  common  is  described. 

The  most  important  common  is  the  labeled  common  MODEL  which 
contains  all  input  and  output  of  the  alternate  energy  models.  The  common 
block  structure  is  illustrated  in  Figure  4-4.  Essentially,  each  model  has 
a  block  of  storage  equivalenced  to  the  AMODEL  array  in  labeled  common 
MODEL.  This  storage  is  positional  for  each  alternate  energy  model 
according  to  the  MTYPE  function. 

Most  of  the  communication  between  subroutines  is  done  through 
labeled  common  blocks.  A  general  description  of  the  use  and  contents  of 
these  common  blocks  is: 

Common/CMMRCL/:  working  storage  for  commercial  energy  parameters. 

These  parameters  are  used  to  store  the  cost  and  amount 
of  commercial  energy  purchased. 

Common/CNTRLS/:  model  and  demand  names,  working  storage  for  number  and 

names  of  models,  and  demands  considered.  This  common 
block  also  contains  values  describing  the  demands  that 
are  satisfied  by  each  alternate  energy  model. 
Common/CQSTS/:  working  storage  for  costs  of  energy  produced  by 

alternate  energy  models  and  commercial  energy  costs. 
Common/ENRGY/:  storage  for  energy  produced  by  alternate  energy 

models.  Each  model  stores  the  energy  produced  hourly 
into  any  of  three  demand  sectors. 

Common/FUELSC/ :  storage  for  site  input  related  to  commercial  energy 

purchases. 
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Common/HEADNG/: 

Common/10/: 

Common /LIN COM/: 

Conrrno  n/MODEL/ : 

Common/NDAYS/: 

Common/SITE/: 

Common/USER/: 


expanded  names  for  demands,  used  in  output  headings, 
variables  related  to  input/output  and  logical  variables 
related  to  end  and  error  conditions, 
input  data  defined  in  INPUT  that  is  transmitted  to  the 
optimization  routines. 

a  large  single  dimension  array  (AMODL'L)  equi valenced  to 
individual  model  parameters  (Figure  4-4). 
storage  for  data  such  as  number  of  days  in  each  month  and 
total  number  of  days  in  the  year. 

storage  for  site  input  data  and  working  demand  arrays, 
storage  for  equality  and  inequality  constraints  and  partials 
of  constraints. 
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SECTION  5 


SUBROUTINE  DESCRIPTION! 


This  section  describes  the  routines  developed  by  Acurex  to 
supplement  the  optimizing  routines  supplied  by  the  Civil  Engineering 
Laboratory. 

Block  Data  BLKDAT :  various  data  elements  used  by  the  code.  Names  of 

models  and  demand  can  be  found  here  in  arrays  MODLS 
and  DMANDS,  respectively. 

Subroutine  BLKOUT:  used  to  print  the  various  yearly  demand, 

insolation,  and  wind  profiles  used  by  the  code. 

Subroutine  COSTCE:  computes  annual  cost  of  commercial  energy  to  meet 

demands  left  by  alternate  energy  models. 

Subroutine  CSUM:  loads  the  BSUM  array  with  values  of  the  equality 

constraints,  and  loads  the  BSUMNT  array  with  values 
of  the  inequality  constraints.  The  code  currently 
has  three  inequality  constraints  on  area,  coal,  and 
refuse  availability.  The  sum  of  area  used  by  the 
alternate  energy  models  must  be  less  or  equal  to 
the  area  available  at  the  site.  The  total  amount 
of  refuse  used  by  the  RDF  models  must  not  exceed 
the  amount  available  at  the  site.  The  amount  of 
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coal  used  by  alternate  energy  models  must  not 
exceed  the  coal  available  at  the  site. 

Subroutine  DETEXC:  determines  excess  energy  produced  by  each  model. 

The  models  are  ordered  according  to  energy  cost  so 
that  the  cheapest  energy  is  used  first. 

Main  Program  ENSITE:  driver  routine  for  the  Navy  Energy  Siting  Code. 

Subroutine  ERROUT:  contains  various  input  error  messages  used  by 

code.  Control  is  not  returned  to  calling  routine 
and  code  terminates. 

Subroutine  EVAL;  executive  driver  for  the  individual  alternate 

energy  models.  EVAL  calls  LOADDM  to  load  demands 
into  working  array  and  then  calls  each  model  to 
determine  energy  produced  and  compute  costs. 

Subroutine  EVALn:  corresponds  to  the  various  alternate  energy  models 

called  by  the  EVAL  routine.  Each  EVAL  ,  where  n 

n 

corresponds  to  the  model  numbers  (see  alternate 
energy  models  list  Table  1-1),  models  a  particular 
type  of  alternate  energy  system.  Given  the  system 
size  (e.g.,  ton/day,  10,000  ft  solar  panels), 
the  models  determine  the  amount  of  energy  delivered 
to  a  particular  demand  or  set  of  demands 
(cogeneration).  Also,  capital  cost,  maintenance 
cost,  fuel  cost,  and  area  required  are  computed. 

Subroutine  FCHART:  uses  f-chart  method  to  determine  the  amount  of 

energy  supplied  by  the  solar  thermal  model. 
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Subroutine  INMODL: 


Subroutine  INPTCL: 


Subroutine  INPT08: 


Subroutine  INPT09: 


Subroutine  INPTXX: 


Subrouti ne  INPUT: 


XL ( 50 )  - 
X8 ( 50 )  - 
XH( 50)  - 


recognizes  model  name  in  the  input  deck  and  calls 
MOOLIN  to  input  model  data. 

a  special  input  routine  called  in  addition  to  INPUT 
(see  description)  for  all  models  which  use  coal  as 
a  fuel.  It  sets  the  capital  cost  factor  and  fuel 
cost  to  correct  values  depending  on  coal  quality 
(high  or  low  sulfur  content). 

a  special  input  routine  called  in  addition  to  INPUT 
(see  description)  when  geothermal  model  EVAL08  is  a 
candidate  model.  It  computes  maximum  potential 
thermal  energy. 

a  special  input  routine  called,  in  addition  to 
INPUT  (see  description),  when  geothermal  model 
EVAL09  is  a  candidate  model.  It  computes  maximum 
potential  electrical  energy, 
called  by  INPUT  to  load  minimum,  maximum  and 
starting  X  values  (number  of  systems)  for  the 
alternate  energy  models. 

initializes  data  for  starting  the  optimization 
process.  It  loads  this  data  into  the  common  block, 
LINCOM,  which  transmits  the  data  to  the 
optimization  routines.  The  variables  used  and 
their  meaning  are  as  follows: 
lower  bound  estimate  for  the  number  of  systems 
starting  values  for  the  number  of  systems 
upper  bound  estimate  for  the  number  of  systems 
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r 


KNT(50)  - 

KON  ( 50 )  - 
N 

CMM 

NLINEQ  - 
NNOTEQ  - 
TOLCON  - 

I  OUT 

ITER 

ITERMX  - 


Subrouti ne  INS  I TE : 


Subroutine  ITRPRT: 


array  used  by  the  system  only 

array  used  by  the  system  only 

number  of  models  being  considered 

-1  minimize  cost 

number  of  equality  constraints 

number  of  inequality  constraints 

allowable  violation  of  constraints,  approximately 

( ( XH ( 1 ) )  -  XL(1))  *  10'7,  never  0.0 

output  flag:  -1  --  answer  only,  0,  1  or  2  --  more 

debugging  information 

used  by  the  system  only 

used  by  the  system  only 

It  also  calls  INPTXX  to  load  the  starting  values 
XL,  XH,  XB  for  each  alternate  energy  model.  If  a 
particular  model,  for  example  geothermal,  requires 
additional  calculations  at  this  stage,  INPUT  calls 
additional  input  subroutines  designed  specifically 
for  that  model . 

reads  the  various  site  dependent  data.  That  data 
includes  the  energy  demand  profiles,  weather 
profiles,  commercial  energy  data  set,  and  general 
data  set. 

used  to  print  progress  of  optimization  during  the 
iterative  process.  This  print  output  is  initiated 
when  IOUT  is  greater  than  -1. 
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Subroutine  LOADDM: 


Function  MATBIT: 


Subroutine  MATCH: 


Subroutine  MAXVAL: 


Subroutine  MODLIN: 


Subroutine  MODOUT: 


Function  MTYPE: 


called  by  EVAL  to  load  energy  demand  profiles  into 
the  working  demand  arrays. 

called  to  determine  the  demand  a  particular  model 
satisfies.  Each  model  has  a  value  in  the  ISTSFY 
array  which,  when  broken  down  into  bit 
representation,  will  identify  the  demand(s)  the 
model  will  satisfy. 

called  by  RUNSPC  to  determine  from  input  which 

demand  or  model  is  being  requested. 

a  general  utility  routine  that  determines  the 

maximum  value  contained  i.:  an  array. 

called  by  INMODL  to  read  data  for  the  model  named 

on  lead  card  of  model  input  data  set. 

called  to  output  the  model  data  in  a  readable 

format . 

used  to  determine  the  subscript  of  model  data  when 
given  the  model  number  (see  Table  2-1  for  model 
numbers).  This  is  done  because  of  the  complex 
structure  of  the  model  data  contained  in  the 
labeled  common  MODEL.  The  user  should  refer  to 
Section  4.4  (common  block  storage)  and  Figure  4-4 
(common  block  structure  for  alternate  energy  model 
parameters). 


Subroutine  OPTMIZ 
(NONLIN): 


the  gradient  projection  optimization  technique, 
supplied  by  the  Civil  Engineering  Laboratory. 


Subroutine  ORDER: 

Subroutine  PARTL(X): 

Subroutine  PRTLIN: 

Subroutine  PRTSUM: 

Subroutine  RAPUP: 

Subroutine  RUNSPC: 

Subroutine  SCRIPT; 

Subroutine  SETOUT: 
Function  UAC: 

Subroutine  UPDATE: 


used  to  reorder  model  calling  sequence  according 
to  cost. 

loads  the  FPC  array  with  the  partial  derivatives 
of  the  constraints,  with  respect  to  the  number  of 
systems.  For  example,  FPC ( I , J )  =  partial 
derivative  of  constraint  J  with  respect  to  number 
of  systems  I. 

called  by  MODOUT  to  print  various  lines  of 
alternate  energy  model  output, 
called  to  give  a  sunmary  printout  after  the 
optimum  is  reached. 

duplicates  last  call  to  EVAL  with  models 
reordered. 

reads  run  specification  input  data  deck  and 
determines  whether  data  are  available  for  models 
requested. 

returns  the  proper  subscripts  to  be  used  in 

referencing  working  arrays  when  given  the  month, 

hour,  model,  and  demand. 

called  to  set  up  output  for  ITRPRT. 

computes  uniform  annual  cost,  given  model  initial 

capital  costs,  and  inflation  rates. 

sets  up  the  last  call  to  EVAL,  applying  the 

maximum  number  of  systems  for  each  demand  to  the 

other  models  in  that  demand.  This  is  done  so 

that  cost  comparisons  can  be  made  across  models. 

An  additional  summary  printout  shows  the  results 

of  UPDATE. 
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SECTION  6 


SAMPLE  RUN 


Presented  in  this  section  is  a  sample  run  which  was  run  on  a  CDC 
7600  Computer.  For  this  run  the  following  is  presented: 

•  A  brief  description  of  the  nature  of  the  run  and  solution 

•  A  listing  of  the  input  data  deck 

0  A  listing  of  the  output  generated 

The  goal  of  this  run  was  to  determine  the  optimum  mix  of  alternate  energy 
systems  and  commercially  produced  energy  at  the  Norfolk  Naval  Base  which 
satisfies  the  three  energy  demands:  space  heating  and  hot  water  (SPCHTG), 
process  steam  (PROSTM),  and  electricity  (ELECTR). 

The  candidate  models  for  this  run  are  refuse  derived  fuel  for  each 
demand  (RDFHTG,  RDFSTM,  RDFELE),  fluidized  bed  combustion  for  each  demand 
(FBCHTG,  FBCSTM,  FBCELE) ;  conventional  coal  for  each  demand  (CCLHTG, 
CCLSTM,  CCLELE);  plus  the  co-generation  model  (CCLCOG)  for  electricity  and 
steam,  solar  thermal  for  heating  (SLTHTG),  photovoltaic  for  electricity 
(PHVEIE),  and  three  wind  electric  models  sized  at  5  kW,  200  kW,  1500  kW, 
(WD5,  WD200,  WD1500) . 

The  input  deck  and  results  are  listed.  Note  that  the  print  flag 
for  this  run  was  set  at  -1  (card  320  of  input  deck). 

The  output  generated  shows  all  demand,  insolation  and  wind  profiles 
input.  Next,  the  geothermal,  coirmercial ,  and  general  site  data  are 
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printed.  Once  all  site  data  is  printed,  each  of  the  model's  data  is 
printed  as  it  is  read  from  the  input  deck.  Then  the  run  specification  is 
read,  checked,  and  printed.  The  demand  numbers  and  model  numbers  are 
printed  for  those  demands  and  models  being  considered.  Each  of  the 
model's  lower  bound,  starting  point  and  upper  bound  are  listed  in  order  of 
the  models  input  sequence.  The  optimization  process  is  not  printed  when  a 
print  option  flag  of  -1  is  used.  The  end  of  the  optimization  is  indicated 
by  the  printing  of  "BEST  RESULTS  OBTAINED,"  a  summary  printout  follows. 

The  summary  printout  shows  the  demand  left  to  be  satisfied  by 
oil-fired  boilers  for  the  heating  and  process  steam  demand  and  purchased 
electricity.  Then  the  energy  produced  by  each  model  for  the  optimum  mix 
and  a  summary  printout  for  each  demand,  indicating  the  models  selected 
(optimum  mix)  is  listed.  Finally,  for  cost  comparison,  a  summary  printout 
for  each  demand  is  listed  using  the  same  number  of  systems  for  each  model 
within  a  demand  sector. 

The  results  for  the  sample  run  show  the  optimum  mix  to  be: 

•  Space  heating  and  hot  water 

—  Fluidized  bed  combustion:  295.4  (tons/day) 

•  Process  steam 

--  Refuse  derived  fuel:  6.7  (tons/day) 

—  Fluidized  bed  combustion:  70.95  (tons/day) 

--  Conventional  ccal  cogeneration:  206.9  (tons/day) 

Electricity 

—  Fluidized  bed  combustion:  553.4  (tons/day) 

—  Conventional  coal  cogeneration:  206.9  (tons/day) 
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NOTE:  Cogeneration  competes  in  two  demand  sectors.  Thus,  206.9 
(tons/day)  is  the  total  amount  of  coal  consumed:  206.9  (tons/day)  is  not 
consumed  in  each  demand  sector. 


i 
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NOTE:  Cogeneration  competes  in  two  demand  sectors.  Thus,  206.9 
(tons/day)  is  the  total  amount  of  coal  consumed:  206.9  (tons/day)  is  not 
consumed  in  each  demand  sector. 
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P.£H&T  VHuOUCCO  BY  ALL  MODELS  TO  SATISFY  -  SPACE  HEATING  AND  MOT  WATER 

AN  •  BEFORE  NUlibEH  INDICATES  AN  EXCESS  OF 
tNCRGT  WAS  PRODUCED  BY  THE  MODEL 
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f-iiLOwING  SUMMARY  USES  THE  SAME  NUMBER  OF  SYSTEMS  IN  EACH  DEMAND  SECTOR 
fjR  COST  COMPARISONS.  NOTE  THAT  COGENERATION  MODELS  MAY  BE  THE  SAME  IN  ONLY 
ONE  DEMAND  SECTOR. 


ENERGY  SOURCE  KUMBCR  OF  DELIVERED  INITIAL  ANNUALIZED  DELIVERED  AREA  REQUIRED  EXCESS  PRODUCED 

STSTE"  UMTS  SYSTEMS  ENERGY  CAPITAL  COST  COST  ENERGY  COST  ENEKGY  ENERGY  COST 

msTui  (n*>  (MU  is/mbtui  ift**2»  imoIui  <*/mbtui 


6-36 


